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Abstract. The Jordan-Wigner transformation is applied to study the ground state properties and dimeriza- 
tion transition in the Ji — J2 XXZ chain. We consider different solutions of the mean-field approximation 
for the transformed Hamiltonian. Ground state energy and the static structure factor are compared with 
complementary exact diagonalization and good agreement is found near the limit of the Majumdar-Ghosh 
model. Furthermore, the ground state phase diagram is discussed within the mean-field theory. In particu- 
lar, we show that an incommensurate ground state is absent for large J2 in a fully self-consistent mean-field 
analysis. 
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1 Introduction 

Low-dimensional quantum spin systems with competing 
interaction are a field of diverse theoretical studies pQ. 
Recent progress in material synthesis has brought up sev- 
eral transition metal oxides realizing such spin systems. 

1 In particular, the s = 1/2 spin chain with frustrated next- 
nearest neighbor interaction has been found to model the 
quasi one-dimensional copper oxide material CuGeC>3 |2I3| . 
The J\ — J2 Heisenberg chain is well known to display a 
quantum phase transition from a gapless, translationally 
invariant ground state with algebraic spin-correlations to 

I a dimer state H] at J 2 ~ 0.24Ji |5I6| . At the Majumdar- 
Ghosh point |2j, i.e. at J 2 = Ji/2, the ground state is a 

' doubly degenerate dimer product of singlet pairs on neigh- 
boring sites. In the vicinity of the Majumdar-Ghosh point 
there is a transition [H] to an incommensurate phase for 
larger J 2 |5I!)| . Low-energy field theories |5I10) are not suit- 
able for describing this transition since it occurs at a very 
short correlation length of about one lattice spacing. Also 
semiclassical approaches |11I12) are not applicable to this 
commensurate-incommensurate transition in the J\ — J 2 
chain with s = 1/2, since they fail to reproduce e.g. the 
solitonic nature of the fundamental excitations |13| . 

The Jordan-Wigner transformation is another tech- 
nique which allows for an approximate analytic analysis of 
low-dimensional quantum spin systems (see e.g. Ref. [jl] 
for a recent review). It has also been applied to the J\ — J2 
chain |15I16| . however only for a restricted set of mean- 
field configurations. While Ref. ^BJ reported incommensu- 
rability for larger J 2 , this work did not consider dimerized 
states and thus fails to reproduce the gap which is known 



to be relevant in this region (see e.g. [H])- This motivates 
us to perform a systematic study of the spin- 1/2 J± — J2 
chain using the Jordan-Wigner transformation and con- 
sidering all possible mean-field configurations. Since the 
Jordan-Wigner transformation is known to be exact for 
the nearest- neighbor XY chain | 17) . we will further con- 
sider the case of a general XXZ anisotropy. Note that the 
effects of the anisotropy in the spin- 1/2 J\~ J2 chain have 
been studied systematically either numerically |18I19| or 
using field theoretical approaches |4I2UI21| , but not in pre- 
vious Jordan-Wigner treatments |15I16| . 

This paper is organized as follows: In section[5]we per- 
form a mean-field analysis of the Jordan-Wigner trans- 
formed 22 23 Hamiltonian. Remarkably, the Majumdar- 
Ghosh state |7I24) is recovered exactly by our approach. 
Section [3] presents results obtained from a numerical so- 
lution of the mean-field equations and a comparison with 
complementary exact diagonalization data. Finally, sec- 
tion ^summarizes our results concerning in particular the 
applicability of the Jordan-Wigner approach to the differ- 
ent parameter regimes. 



2 Mean-field analysis 

In the following we consider the frustrated spin-1/2 XXZ 
chain of L spins (L — > 00): 



1=1 



+J 2 (sfsf +2 + sfsf +2 + Asfsf +2 



(1) 
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Here J\ (J 2 ) > is the nearest (next-nearest) neighbor 
coupling, A is the exchange anisotropy and sf are spin-1/2 
operators. By the Jordan- Wigner transformation 22 23 
the spin Hamiltonian Q is mapped onto a model of in- 
teracting spinless fermions c 



(+). 



1=1 



h.c. ) + AJi [ cfci 



c[ +1 ci +1 - 



AJ 2 c+ci 



1 



C l+2 Cl+2~ 



(2) 



The first term corresponds to the XY nearest-neighbor in- 
teraction and can be treated rigorously 1251 . The remain- 
ing terms are the z component of the nearest neighbor and 
the complete next-nearest neighbor exchange. They rep- 
resent four-fermion interactions which render the model 
non-integrable. We will therefore resort to a mean-field ap- 
proximation preserving all pair correlations of type (c+c m ) 
in the factorization of the interaction terms: 



1 



= \ s i ) = K c i C V - 



fji 

Ai = (c+q + i) - 
A = (c/"q +2 ) = (c? +2 ci). 



(3) 



These contractions are related to single-site, nearest neigh- 
bor, and next-nearest neighbor spin-spin correlation func- 
tions. The self-consistent determination of these contrac- 
tions follows Ref. |15|. where however only the case of 
gi = Di = had been considered. We find the following 
phases: 

i) paramagnetic {homogeneous) 26 : A\ = A = — — , gi — 
A = 

ii) uniform antiferromagnetic (AFM+uni ): gi = (—l) l g, 
Ai = A, and Di = D 

iii) staggered antiferromagnetic (AFM+stag.): in contrast 
to ii) Di is staggered, i.e. Di = ( — l) l D 

iv) alternating nearest neighbor hopping (dimer): Ai = 
A + (-l) l S, 

The self-consistency equations for phase iv) are discussed 
below. Here, we should only note that in all cases we find 
no uniform contribution to the next-nearest neighbor cor- 
relation Di. A staggered next-nearest neighbor contrac- 
tion is induced by the antiferromagnetic order and tends 
to zero when g — 0. Therefore, the current mean-field 
approach is not suited to treat the limit of two weakly 
coupled chains (J\ — ► 0) where the correlation between 
the next-nearest neighbor spins becomes strongest. 

For case iv) , and after Fourier transformation the mean- 
field Hamiltonian reads 



H = 



- 7T < k < 7T 



where e k = (Ji - 2A(AJ X - 2J 2 )) cos(fc), f k = 25(AJ 1 + 
2J 2 ) sin (A:). The Hamiltonian can be diagonalized by the 



unitary transformation 



Cfe 



Ck±TT 



1 + 1 

1 - i 



(cos(a{/2) r) k - sin(a{/2) r/ k±n ) , 
(sin(a{/2) rj k + cos(a{/2) r)k±n) , 



where cos a{ = e k j \J ef+/7. This leads to the free fermion 
model 



(5) 



:<k<i 



with the spectrum: \ k = \J z\ + f k ■ For the selfconsistent 
equations for A, 5 we find: 



A = ~ T, 



5 = T, 



dk cos(fc) 



dk sin(fc) 



V4 

fk 



Pk 



(G) 



The mean-field approximation to the ground-state energy 
is obtained from the ground state expectation value of the 
fermion Hamiltonian including contractions only up 
to quadratic order. For the dimer phase iv) the ground 
state energy per site is 



e=J x A- (AJi - 2J 2 )A 2 - {AJi + 2 J 2 )S 2 



(7) 



We note that the dimer order parameter [5] <i=(S2i-iS2 S :) — 
(S2iS2i+i) corresponds to A(A— 1)5 within our treatment. 

So far, we have assumed that the mean fields are real. 
However, it is known that an easy-plane anisotropy A < 1 
induces a chiral phase for large J 2 |1UI27| (note that the ex- 
istence of such a chiral phase had been verified numerically 
in a different ladder model j^Hl prior to [23). The chiral 
order parameter Ki — sfs^ +1 — s^sf +1 corresponds to the 
imaginary part of the fermion nearest-neighbor contrac- 
tion 9rm4/. Therefore, the elementary contraction should 
be allowed to be complex. In this case the ground state 
energy reads: 



,hA R - (AJ! - 2J 2 )A\ - {A.h + 2J 2 )5 2 R 
-(AJi + 2J 2 )Aj - {AJi - 2J 2 )5\ . 



(8) 



Here Ar(J) is the real (imaginary) part of A, and <5r(j) 
is the real (imaginary) part of 5. It is evident that for 
large J 2 / J\ a non-zero Aj lowers the ground state energy 
in contrast to 5[. Hence, we have considered mean- field 
solutions with non-zero Ar, Ai, and 5r. However, numer- 
ical solution of the mean-field equations did not yield any 
complex solution. 

The Majumdar-Ghosh point [7I23| J 2 = Ji/2 per- 
mits to check the consistency of the mean-field approx- 
imation. At this point the ground state is doubly degen- 
erate and consists of singlet pairs on neighboring sites 
njla^ + MZ], where [21^1,21] denotes one singlet. This 
can be related to the fermion language by recalling that 
S z = +(— )l/2 on site I corresponds to a filled (empty) 
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Fig. 1. Ground-state energy per site as a function of J2/J1 
for the Ji — J2 XV chain (upper panel) and the Heisenberg 
chain (lower panel): the circles are exact diagonalization data 
extrapolated to the infinite system; lines correspond to the 
different mean-field solutions for the ground state. 



site /. Therefore (c^ Tl — c^lO) creates a singlet bond 
on neighboring sites 2Z =F 1, 21, and the Majumdar- Ghosh 

states can be represented [25] as ( c 2i=ci — c ti We 

will now argue that this exact state is also obtained from 
the mean- field solution. For J 2 = J\/2 configuration iv) 
and Eqn. © yield A = —1/4, 6 = ±1/4. Inserting this 
solution into the mean-field Hamiltonian one gets 



H = 



1=1 



A) 



where ^ = -^{c 2 i T i+c 2 i), r)\ s) = -^(c 2 i T i-c 2 i) for 6 = 



v'2 



±|. ??j and jj| 8 ' (i7j and rj; ) create and annihilate 
a singlet (triplet) on nearest neighbors. This shows that 
the mean-field solution is the Majumdar-Ghosh state in 

fermionic representation Ilila 2 ? 7/ S ^ + |0) with the ground 
state energy per spin e = — J%{2 + A)/8. 



3 Results 

Fig. n shows results of a numerical solution of the self- 
consistency equations for configurations i)-iv). It can be 
seen that all solutions exist for positive J 2 . To check the 
accuracy of these results, we have performed exact diag- 
onalization of the Hamiltonian for L < 32 sites (since 
this type of computations has a long tradition |30j . sim- 
ilar results can also be found in the literature, in partic- 
ular for A = 1). We have performed a finite-size extrap- 
olation using the Vanden-Broeck-Schwartz-algorithm [3"T| 




Fig. 2. Phase diagram for the J1 — J2 XX Z chain as a function 
of the anisotropy A: the solid line is the result of the mean-field 
approximation, the circles are results of exact diagonalization 



with q;vbs — — 1 and, for most cases, L = 12, 16, . . . 32. 
The extrapolated results for the ground-state energy per 
site e are shown by the circles in Fig. ^ Errors are esti- 
mated to be smaller than the size of the symbols. For J 2 = 
0, our extrapolated values agree to better than 10 _6 Ji 
with the exact answer e/Ji = —1/tt » —0.318310 and 
1/4 - ln(2) w -0.443147 for A = and 1, respectively. 

Apart from the exact ground state at the Majumdar- 
Ghosh point we may contrast the mean-field solutions 
against other known results. We start from the XY nearest- 
neighbor model (A = 0, J 2 = 0) where the Jordan- Wigner 
transformation yields the exact answer JJj- For small J 2 j J\ 
the results of all approximations are very close (upper 
panel of Fig.^) and all order parameters grow very slowly. 
In fact, a detailed analysis of © for A, S shows that 
J 2 > J2,dim = —AJi/2 is the condition of a non-zero solu- 
tion for S. The approximate solution for small J 2 — J2,dim 
above this point gives 



e( J 2 — J2,dim) 



exp 



4( J 2 — J2,dim) 



(10) 



where a — J\ — ^ ( J 2 + J2.dim) • The excitation gap appears 
at k = 7r/2 and increases as 4(J2 — J2,dim)<5. Note that 
close to the critical point J 2c bosonization yields a very 

similar form j 1I5| 5 ~ exp p~r^j~j ) an d a gap which is 

proportional to S 2 . 

For an analysis of the phase diagram we compare the 
ground state energies of all phases. In the regime of small 
frustration, the " AFM '+ stag, ." solution has the lowest en- 
ergy and at some point a c = J 2c / J\ crosses with the so- 
lution for the dimer state (see Fig. Q] for A — 1). This 
can be identified with a first-order phase transition. Note 
that at A = 1 the "AFM+uni." and the dimer solutions 
cross very close to an early estimate for the critical point 
J 2 = J1/6. Taking into account the "AFM+stag." so- 
lution shifts this crossing point, thus yielding very good 
agreement with numerical results |] J2 C ~ 0.242Ji for 
A = 1. Fig. |21 shows this crossing point as a function of 
the anisotropy A. Although the critical value of J 2 is very 
close to numerical results for A w 1, the position of the 
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Fig. 3. Excitation gap in the mean-field Hamiltonian as a 
function of J2/J1 for the Ji — Ji XY chain (upper panel) and 
the Heisenberg chain (lower panel). 




Fig. 4. The structure factor S zz (q) as a function of q for the 
Ji — Ji XY chain (upper panel) and Heisenberg chain (lower 
panel) at J% = 5Ji/ll. The line is the result of the dimer 
mean-field solution; symbols are exact diagonalization data for 
different system sizes from 20 to 32 sites. 



crossing decreases with decreasing A whereas the exact 
diagonalization data [18119] exhibit an increase of Ji c . 

The excitation gaps of the different mean-field solu- 
tions are depicted in Fig. The "uniform!'' solution is 
always gapless and is not displayed. The "dimer" and 
"AFM+stag." solution have a maximum of the gap near 
the Majumdar- Ghosh point which is consistent with nu- 
merical results |5]. The position of the minimum in the 
excitation spectrum of the "dimer" solution jumps from 
k = ±7r/2 to k = at the Majumdar- Ghosh point. The 
minimum of the u AFM+stag." solution shifts continuously 
for J 2 > Ji/2. 

The 2z-structure factor S zz (q) — Ylf=i ex P(*90( s n s n+z) 
can be expressed through the four-fermion correlation func- 
tion. For the dimer state we get the following result in the 
mean-field approximation: 

S zz (q) = --— f dk e ^k+ q + fkf k+q (n) 

In the limit of the Majumdar-Ghosh model the exact re- 
sult [H] S zz (q) = ±(l-cos(q)) is recovered. The result ((TTJ 
for the zz static structure factor at J2/J1 = 5/11 is shown 
in Fig.0]as a full line in comparison with numerical results 
(open symbols). Both approaches yield a broad maximum 
at the boundary of the Brillouin zone q = ±7r, signify- 
ing antiferromagnetic correlations with a short correlation 
length |5I!J |. Small deviations are observed at q = ±7r, 
where the numerical data exhibits also the biggest finite- 
size effects. 



4 Conclusions 

To summarize we have applied the Jordan- Wigner trans- 
formation to the spin-1/2 J\ — J 2 XX Z chain. We have 
performed a complete analysis of the possible mean-field 
states and found that a dimerized state has lowest en- 
ergy for J2 larger than some critical value. The ground- 
state energies are in good agreement with exact diagonal- 
ization data and we even recover the exact result at the 
Majumdar-Ghosh point [7]. 

The location of the critical point is obtained from the 
crossing between the energies of the antiferromagnetically 
ordered state iii) and the dimer state iv). This yields good 
agreement with numerical results for 4»1, but does not 
follow the anisotropy dependence for small A correctly. 
This discrepancy may be related to the following facts. 
Firstly, spin-spin correlations should decay as a power-law 
below Ji c @|, while the mean-field theory treats this phase 
as antiferromagnetically ordered. Secondly, in the mean- 
field scenario the transition is generically of first order, 
not of the established infinite order. Just at A — 0, we 
find Ji c = J2,dim = and recover the infinite order phase 
transition to the dimer state, see l|10|) . 

Furthermore, we have calculated the zz static struc- 
ture factor for the "dimer'' mean-field approximation. The 
result coincides with the exact one for the Majumdar- 
Ghosh point and satisfactory agreement is found with nu- 
merical results in the dimerized phase below the Majumdar- 
Ghosh point. 

It should be mentioned that there are problems with 
the present approach for large In particular, the mean- 
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field dimer state does not reveal any incommensurability, 
as signaled e.g. by a shift of the maxima of the static struc- 
ture factor beyond the Majumdar- Ghosh point This is 
an important difference with Ref. which obtained in- 
commensurability from a mean-field approach. However, 
the solution of Ref. fTf)] does not correspond to the ground 
state. Remarkably, the absence of incommensurability is 
a special property of the dimerized mean-field solution 
which leads to a cancellation of the next-nearest neighbor 
hopping processes in the fermionic picture. Consequently, 
alternative approaches are needed to describe the region 
J 2 > Ji/2. 
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